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Abstract 

We study the dynamics of spontaneous generation of coherence and photon spin-qubit entangle- 
ment in a A system with non-degenerate lower levels. The cases of entanglement in frequency only 
and frequency and polarization are compared and the reduced density matrix and entanglement 
entropy are analyzed. We explore in detail how which-path information manifest when the energy 
difference between the qubit states is larger than the linewidth of the excited state suppresses co- 
herence. A framework is provided to describe the dynamics of spontaneous generation of coherence 
and (ideal) photodetection obtaining the post measurement qubit density matrix. A simple model 
of photodetection with a quantum eraser to suppress which-path information in the detection mea- 
surement is implemented. It is found that such quantum eraser purifies the qubit density matrix 
after photodetection, our results are in agreement with those reported in recent experiments. 
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I. INTRODUCTION 



Quantum entanglement has evolved from being a paradoxical aspect of quantum 



mechanics [1| to becoming a resource for quantum computing and r 
with potential for technological breakthroughs in these areas)* 
ments demonstrated photon entanglement with single atoms [8|- 



luantum information [2 
Several recent experi- 



atomic ensembles 11 



optical cavities [151. 116 



long-distance entanglement between qubits|12Hl4||. and tunable ion^photon entanglement in 



Along with atom-photon entanglement 0, [8 10 1 , and entanglement 
in cavity quantum electrodynamicsl [l~7t] recent proposals suggested electron spin-photon en- 
tanglement in quantum dots as platforms for entanglement between distant spins[l8|]. Spin- 
photon entanglement could be the pathway towards implementation of quantum networks 
among distant nodes jl, 0, 3]. Remarkable experiments demonstrated the realization of 
entanglement between the polarization of a single optical photon and an electronic spin 
qubit in nitrogen vacancy (NV) centers in diamond 19j and more recently the demonstra- 
tion of entanglement between a single electron spin and a photon in a quantum dot has 
been reported[20-22||. A main paradigm in many of these experiments is that of sponta- 
neous generation of coherence 23M25| in a type-II or A system, namely a situation in which 
spontaneous emission from a single excited state via a two-channel decay to degenerate or 
non-degenerate lower levels results in coherence between these two states. Spontaneous gen- 
eration of electron spin coherence has also been observed from the radiative decay of charged 



excitons (trions) in quantum dots [26 



These experimental efforts are paving the way towards the implementation of atom- 
photon or spin-photon entanglement as potential platforms for quantum information and 
quantum computing protocols and networks [3_ll4], |27| , motivating a theoretical effort seeking 
a deeper understanding of these processes [24. |28|-|30|]. 

Although there have been some recent studies of the dynamics of spontaneously generated 



coherence [24l. l25l. |30| many important aspects merit further investigation. 

Our main goal in this article is to provide a more complete theoretical study of the exper- 
imental results reported in ref. 19 [but that apply more generally to current experiments on 
spin-qubit-photon entanglement 20-22|| from spontaneous generation of coherence as men- 
tioned above. With this aim, we focus on the following aspects: 1) to provide a treatment of 
the dynamics of spontaneous generation of coherence, entanglement both in frequency and 
polarization and photodetection within a single framework consistent with causality (3l|. 
2) to study the entanglement entropy of reduced spin-qubit density matrices after tracing 
over the radiation degrees of freedom for photon-qubit entanglement both in frequency and 
polarization, of particular interest when spontaneous emission produces polarized photons 
which are measured by projection on differen polarization states 3) to analyze in detail how 
which- path information affects coherence, in particular within the setting of the experi- 
ment in ref. [l9|, predicting the time dependence of conditional probabilities when which 
path information is present. 4) To implement a model for a "quantum eraser" jH, 33 1 within 
the framework of photodetection a Id Glauber [34 -36] so as to erase which path information 
in the photodetection process. An important result of this treatment is that "quantum 
erasing" "which- path" information leads to the purification of the qubit state confirming 
the experimental results of ref. 19j, |20| and bolstering the arguments on "quantum erasing" 
in these references. We obtain a conditional probability in complete agreement with the 
experimental results of ref. 

[ll- 
Our study differs from and complements recent theoretical treatments of spontaneous 
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generation of coherence 24], |30] in that we analyze both frequency and polarization entan- 
glement, which-path decoherence, the spin-qubit entanglement entropy and incorporate a 
Glauber model of broadband photodetection[34,35,36] in a unified manner with the treat- 
ment of spontaneous emission. This treatment directly builds in causality in the spontaneous 
emission/photodetection process[31], leads to detailed understanding of how which path in- 
formation affects coherence, and allows to model a quantum eraser[32,33] consistently within 
the broadband photodetector model. This approach is different from that advocated in a 
recent article|30( where the photodetector is modeled with a collection of two-state atoms 
spread over some distance where the excited jstate features a short lifetime. Furthermore 
our study also differs from those of refs.[2^ 



30| in that it shows how the implementation 



of a "quantum eraser" leads to the purification of the qubit density matrix upon photode- 
tection and yields a result for the conditional probability in complete agreement with the 
experimental findings in ref. [19]. 



II. DYNAMICS OF ENTANGLEMENT VIA SPONTANEOUS DECAY. 

We consider a A-system with one excited state \A) and two Zeeman split non-degenerate 
lower levels interacting with the electromagnetic field in the dipole and rotating wave approx- 
imations. The degenerate case can be obtained straightforwardly. We refer to the two-lower 
state levels | ± 1) as a spin- qubit. The cases in which there is photon-qubit entanglement 
in frequency only and in frequency and polarization are studied separately and compared. 



A. Entanglement in frequency only: 

We first consider the case when the dipole matrix elements are independent of the polar- 
ization of the photon and for simplicity we only consider one polarization to establish contact 
with the results of ref. [24]]. This case leads to qubit-photon entanglement in frequency only, 
and generalization to two polarizations is straightforward. The total Hamiltonian for the 
three level system is given by 

H = H A + H R + H AR , (II.l) 

where 

H A = E A \A)(A\+E + \ + l)(+l\ + E.\-l)(-l\ ; H R = Y J ^a\a % . {11.2) 

k 

The interaction Hamiltonian in the interaction picture and in the rotating wave approxima- 
tion is given by 



H AR (t) = 1 9k at [| + 1) (A\ e^"^ + | - 1) (A\ e 



i(k-Q-)t 



+ h.c.} (11.3) 



k 

where 



Q ± = E A -E ± ; gk = -ijJ^D, (II.4) 

here V is the volume and D is the dipole matrix element neglecting polarization degrees of 
freedom. 
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Consider that at time t = the initial state is 

|¥(0)> = \A) |0 7 > 



(II.5) 



where |0 7 ) is the radiation vacuum state, and write in the interaction picture the time 
evolved state as 



|*(t)>, = c A (t)\A)\o 7 ) + \W hAt) I + 1) + CkAt) I - 1) 



(II.6) 



The coefficients obey the following equations (in obvious notation) 

C A (t) = -z(A;0 7 |iJ^(t)|%;+l)C , fc , + (t)-2(A;0 7 |F Ai? (t)|lg;-l)C' fc ,_(t) (II.7) 
CkM = -i{l$±l\H AR {t)\A-%)C A {t). (II.8) 

We solve this system of equations with the initial conditions 

C A (0) = 1; (7 fc ,±(0) = 0. (II.9) 



In the Wigner-Weisskopf approximation [34], |37j the coefficients are given by 1 



(11.10) 
(11.11) 

r = r+ + r_ (11.12) 

where the partial widths I~± correspond to the spontaneous decay channels \A; 7 ) — > + 
1); I A) — > |lg)| — 1) respectively, namely 



C A (t) = e"5 
C*,±(*) = Wk 
The level width V is given by 



r± = 2vr ^1(^1^(0)11^; ±1)| 2 S{k - fi±) = 2tt ]T |^| 2 <5(A; - fi ± ) 



2vr 



(11.13) 



In most experimental circumstances, the energy splitting is much smaller than the optical 
frequency of the transitions, namely \Q + — <C Q± in which case it is convenient to write 



Auo 

Q± = Q ± — ; Aw«fl 
and to leading order in Au/Q it follows that 

r + ~ r_ ~ r/2. 



(11.14) 



(11.15) 



In the experiment reported in ref. 19j, it has been verified that the approximation ( 1II.15|) 
is fulfilled in the setting of that experiment. In what follows we will assume that the relation 
(1II.15P holds unless otherwise stated. 



1 We neglect the contribution from the Lamb shift to the energy level E A . 
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We write the spin qubit-photon entangled part of the wavefunction (in the interaction 
picture) (1II.6j) as 

|<M*)> = i [[<ri(t)>| + 1) + \a 2 (t))\ - 1>] (11.16) 
where the single photon wavepackets are given by 

Mt)> = V2j2c k ,+(t)\W ; W 2 (t)) = V2j2c k At)\W (n.i7) 



B. Normalization of photon wavepackets: 



The normalization and orthogonality of the single photon wave-packets is determined by 
the overlaps 

(a a (t)\a b (t)) = 2j2C lb (t)Cl a (t) ; a,b = 1,2. (11.18) 

fc 

Consider the functions 



I — gi(w— f2 a +i~)t 



u — £l a + i ~ 



(11.19) 



in the narrow width limit T fi Q these are sharply localized near u; ~ Q a , straightforward 
contour integration yields 



Ga{w, t)G*Ju, t) duj = 2n 



]_ _ g— i(n a — &p)t g— rt 



(11.20) 



r + z(o a -fy) 

Combining this result with ( 1II.17|II.lip we find consistently with the Wigner-Weisskopf 
approximation 



(cr h 2(t)\a h2 (t)) 



1 



-rt 



(11.21) 



This result, along with the relation between the total and partial decay widths given by 
( 1II.12|) yields the normalization of the \^f ap ) state, 



<tt ap (t)|* flp (t)) 



1-e 



-rt 



(11.22) 



which is a result of unitary time evolution manifest in the Weisskopf-Wigner formulation 
since the total state \^(t))j given by ( \ll.6\i must obey (^(t)\^(t)) = 1. Because |cr 1)2 ) are 
single photon wavepackets, it is straightforward to confirm that the total number of photons 
is given by 

' _ l-e- r *l. (11.23) 



N 1 (t) = (* ap (t) X)4 a * 



Taking r + ~ T_ ~ T/2 under the assumption that Au <C Q, consistent with the ex- 
perimental setup in [19J, it follows that the single photon wavepackets are normalized for 
rt ^> 1 but they are not orthogonal, we find 



M*)ki(*)> 



■y g— iAoit g— rt 



1 + i^ 



(11.24) 
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a result that is in agreement with an observation in ref . [24| for Ft ^> 1 . 

Let us consider the reduced density matrix for the qubit by tracing over the radiation 
field, namely (in the interaction picture) 



p I fo (t)=Tr R \* ap (t)){* ap (jt)\ 
going back to the Schroedinger picture we find 



Pfo(t) 



1-e 



-rt 



+ + e iAaJt 7/(t)| + l)(-l|+h.c. 



(11.25) 



(11.26) 



where 



r,(t) = \ri(t)\e**) 



e -iAu)t e -rt 



(l- c -r.)(l + i^) 



(11.27) 



In the long time limit Tt ^> 1 the coherence is suppressed by the factor 1/y 1 + Auj 2 /T 2 
reflecting the suppression of coherence by "which-path" information. If F ^> Ao; the spectral 
width of the radiation, determined by the lifetime of the excited state, suppresses which path 
information by blurring the energy resolution of the decay channels of the emitted photons 
and coherence is maintained. In the opposite limit Au ^> T the energy difference between 
the lower lying states is resolved and which path information is available in the emission 
spectrum thereby suppressing coherence. This is manifest in the overlap of the photon 
wavepackets (III. 241) in terms of the product of the Lorentzian line shapes for the individual 
channels. 

The main reason for studying the reduced density matrix in the case of frequency en- 
tanglement only is that, as it will be discussed in detail in section (1IIII) photodetection that 
filters horizontal (H) or vertically (V) polarized photons projects the density matrix onto a 
reduced density matrix precisely of the form (1II.26P that contains "which-path" information. 



C. Entanglement in frequency and polarization 

In the experimental situations considered in for atom-photon entanglement and 

in ref. 19) for electron spin-photon entanglement in NV-centers, there are angular momentum 
selection rules in spontaneous decay and the photons emitted are right handed (for \A) — > 
| — 1)) or left handed (for \A) — > | + 1)) circularly polarized as depicted in fig. In 
this case the spin- qubit and the spontaneously emitted photons are entangled both in 
polarization and frequency. Including the polarization of the emitted photons leads to 
several important modifications of the results obtained in the previous case, therefore we 
restore the polarization, momentum and spatial dependence of the dipole matrix elements. 
Although we focus the discussion on the experimental setup of ref. [l9| with NV-centers, the 
results will be more general. 

In this case the total Hamiltonian for the three level A-system interacting with the elec- 
tromagnetic field is given by (III. II) with Ha given in eqn. (III. 21) . but now 

H R = ^4,a g M , (H.28) 

k,\=± 
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FIG. 1: Transitions 

and the interaction Hamiltonian in the interaction picture and in the rotating wave approx- 
imation is given by 



H AR (t) = \ 9k, 
k < 



(x ) a 



.l)( j A| e i(*-n+)* + ^_(f ) a t \-i)(A\e i V°- n -'* t +h.c.\ (11.29) 



where Q± = E A — E±, and 



9U t± \? ) = -i\j — D ± -e % ^ e 



ik-XQ 



(11.30) 



here V is the volume, D± are the dipole matrix elements (±l|d|A) respectively, are 
the left and right handed polarization vectors respectively and x$ is the position of the NV 
center. 

Consider that at time t — the initial state is 



|tf(0)> = \A) |0 7 ) 



(11.31) 



where |0 7 ) is the radiation vacuum state, and following the notation of the previous section 
we write the time evolved state in the interaction picture as 



m))i = c A {t)\A)\%)+Y, 



c t+ {t)\u)\ + i) + c t _{t)\u + )\-i) 



(11.32) 



The coefficients obey the following equations (in obvious notation) 

C A (t) = -i(A;O y \H AR {t)\l n + ;-l)C^_(t) - t(A;0,\H AR (t)\l^; +1) C u (t) (11.33) 
C%Jt) = -i(l^-±l\H AR {t)\A-%)C A {t). (11.34) 

Just as in the previous section we solve this system of equations with the initial conditions 
C A (0) = 1 ; Cr ± (0) = 0, in the Wigner-Weisskopf approximation the coefficients are given 

by 2 



Again we neglect the contribution from the Lamb shift to the energy level E A . 
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e 2 



C A (t) 
C^ ± {t) = ig£ i± {x ) 

The level width T is given by 



i(fe-n ± +ir/2)t 



r\ +r_ 



(11.35) 
(11.36) 

(11.37) 

where the partial widths T + , T_ correspond to the spontaneous decay channels |^4;0 7 ) — > 
— l)j 1^4) |lfc_)l + 1) respectively, namely 

r + = 2vr J2\( A \H AR (o)\h^ +i)l 2 <^ - n+) = 2nJ2 \% + (m 2 S{k - n + ) (I1.38) 

T_ = 2tt ^ 1(^1^(0)11^; +1)| 2 8(k - = 2tt £ |^_(f )| 2 5(A; - . (11.39) 

Just as in the previous case of unpolarized photons, it follows that T± oc |5 ± | 2 n^ but the 
proportionality constants now depend on the angular average of the polarization vectors. 

Now the second term of the wave function ( III. 321) describes an entangled state of circularly 
polarized photons and the spin states of the NV- center, following the literature [g- 10, 2 
we write this second term (in the interaction picture) as 



I'M*)) 



V2 



Mt)> I + 1> + K(t)> I - 1> 



where 



|a T (t)) = v / 2^Q j± (t)|l fe> ) 



(11.40) 
(11.41) 



describe orthogonal circularly polarized single photon wave packets. 

Unlike the results in ref. |24j we do not take the limit Tt ^> 1, in the experimental setting 
of ref. 19[ the lifetime of the excited state is 1/T « 12 ns but the measurements are performed 
during a time interval ~ 10 — 20 ns. 

Borrowing the results from the previous section, we now find 



(a T (t)\a T (t)) 



2r H 



-Ft 



(a + {t)\a-{t)) = 



(11.42) 



where the orthogonality of \a T (t)) is a consequence of the fact that they describe one photon 
wavepackets with orthogonal polarizations. This result, along with the relation between the 
total and partial decay widths given by (111.37j) again yields the normalization of the \^ S p) 
state, 

(11.43) 



which is a result of unitary time evolution and similarly 

N,(t) = (V sp (t)\ a \ x %x M*)) 



-rt 



(11.44) 



fc,A=± 



s 



Just as in the previous section the one-photon wave packets \a±) have unit normalization 
when Ft ^> 1 and T + = T_ = T/2, which is justified when the Zeeman splitting Q + — <C 
and describes the experimental setup of ref.[19j]. The reduced density matrix for 
the spin-qubit can be obtained by tracing over the radiation field just as in the previous 
section flII.25IIII.26p . However, in this case, the orthogonality of the circularly polarized 
wave packets leads to vanishing coherence and a diagonal density matrix that describes a 
statistical mixture given by 



Pfp(t)= r ft\*8 P (t)){*sp(t)\ 



1-e 



-rt 



+ 1)(+1| + |-1)(-1| 



(11.45) 



D. Entanglement entropy: 

As we have seen above spontaneous generation of coherence leads to very different reduced 
density matrices depending on whether photon-qubit entanglement is in frequency and po- 
larization or frequency only. This difference is highlighted by comparing the Von-Neumann 
entanglement entropy in both cases. 

Frequency entanglement only: in this case the total reduced density matrix is 



P(t) 



-rt 



\A)(A\ + Pfo (t) 



(11.46) 



where Pf {t) is given by (1II.26j) which can be diagonalized with the following eigenvectors 
and eigenvalues 



|1> 
|2) 



4=(| + 1) + e~ lv{t) e~ iAuJt \ 
\/2V 1 



1) 



+ l)-e- itp{t) e -' lAult \ - 1) 



1-e 



-rt 



A 2 (t) 



1-e 



-rt 



where r](t) = |^(t)|e 41/3 ^' ) is given by (0X27]), leading to 

p{t) = e- Tt \A){A\ + Ai(f) fx><xj + A 2 (t) |2> <2| . 
The entanglement entropy follows directly, 

S fo {t) = Yte- Tt - A x (t) lnAx(t) - A 2 (t) lnA 2 (t) . 

For rt > 1 



Sfo(<X>) 



1 + 1^0 



111 



1 + k 



Wo 



111 



l+\v(t)\\ (11.47) 
l-\v(t)\] (11-48) 

(11.49) 
(11.50) 

(11.51) 



with 



1+ AuP 



(11.52) 



As Auj/T — > the entanglement entropy vanishes as the asymptotic state is the pure state 
|1) = + | — )j in the opposite limit Au/T ^> 1 where which-path information sup- 

presses coherence it follows that S/ (oo) = ln(2) describing an equal probability statistical 



mixture. 
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Entanglement in frequency and polarization: in this case the total reduced density 
matrix is simply 



\A)(A\ 



1-e 



-rt 



+ !)(+! + | -1)(-1|) 



(11.53) 



as a consequence of the orthogonality of the right and left circular polarized photon 
wavepackets. In this case the entanglement entropy is 



[1 - e~ rt ] In 



1-e 



-rt 



(11.54) 



S fp (t) = Tte- rt 
with the asymptotic value 

S fp {oo) = ln(2) . (11.55) 

The entanglement entropies in both cases are displayed in fig. ([2]) for the parameters of 
the experiment in ref.pl], Au = 2n x 122 MHz ; T = 1/12 ns. 



1.2 




t(ns) 



FIG. 2: Entanglement entropy for the case of entanglement in frequency only Sf a (t) and frequency 
and polarization S fp (t) for Au = 2vr x 122 MHz ; T = l/12ns,r + = T_ = T/2. 



Analytically it can be seen that 



S fp (t) > S fo (t) 



(11.56) 



a relation that is confirmed numerically and confirms the qualitative expectation that the 
entanglement entropy should be larger in the case of entanglement both in frequency and 
polarization. 

The results above were obtained under the assumption that r + = T_ = T/2. If the 
partial widths to the two non-degenerate levels are different the generalized form of the 
entanglement entropy in this case of entanglement in frequency and polarization is given by 



%(t) = ne- rt -^(l-e- r > ^(l-e- rt ) -^(l-e~ r > ^(l-e"") (IL57) 



L r 



where T 
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III. PHOTODETECTION 



We consider a model for a broadband photodetector described by an atom localized at 
position Xd interacting with the radiation field in the dipole approximation a Id Glauber 34- 



361 ] . The Hamiltonian is given by Hp + Hd R where the detector Hamiltonian Hp describes 



a zero energy ground state and a collection of excited states which eventually will be taken 
as a continuum 

H D = v Q \g d ){g d \ + Y,^\4)(ej\ ; u = , (III.l) 

i 

and H^R is the interaction Hamiltonian that describes a dipolar coupling to the radiation 
field with a filter that selects H/V linear polarization states of the radiation field. In the 
rotating wave approximation and in the interaction picture it is given by 

H DR (t) = J2[dj-E P + \x d ;t)\e d )(g d \e^ t + h.c^ ; P = H; V (III.2) 

3 

dj are the dipole matrix elements and 



E p + \x d ;t) = X>V^ ip %p^ d e- m - (HI.3) 



The combined process of spontaneous emission from the NV-center \A) considered to be 
localized at xq = and photodetection by a broadband photodetector localized at x d is now 
described by the total Hamiltonian 

Htot = H A + H R + H AR + H D + H DR , (III.4) 

where the first three terms are given by fill. 11111.3]) . 

Insight into the combined processes and the intermediate states that contribute is gleaned 
in second order in the perturbative expansion with the full interaction Hamiltonian in the 
interaction picture (and in the rotating wave approximation) 

Hj(t) = H AR (t) + H DR (t) (III.5) 

where H AR {t); Hr, R (t) are given by (1II.29j) with x = and (1III.2j) respectively. Consider 
that the initial state is (in obvious notation) 

|*(0)) = L4;0 7 ;<? d ), (III.6) 
in the interaction picture the resulting time dependent state in second order becomes 



!*(*)> 



\ E^dtx + {-if \ \ H^H^dhd^ + 
Jo Jo Jo 



|^(0)). (III. 7) 



To first order only H A r contributes and describes the perturbative spontaneous decay of 
the excited state \A) of the NV-center into the Zeeman split states —1) and |l£_; +1). 
Inserting a complete set of eigenstates of H = H a + Hd + H r it is straightforward to see that 
in the second order contribution the first term Hi fa) describes the spontaneous emission of 
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the circularly polarized photons while the second term Hj(ti) describes the absorption of 
these photons and the photoexcitation of the detector (along with a second order contribu- 
tion from Ha r that yields the original state back). The photodetection probability at time 
t is given by [34 -3^ 

P D (t) = Tr d J2\e*)(e*\p(t), (III.8) 



where the density matrix 

p(f) = |¥(f)><*(f)l, (III.9) 

and the trace in (IIII.8I) is over the detector excited states. 

Our goal is to describe these processes non-perturbatively with a Wigner-Weisskopf de- 
scription that incorporates both processes at once. Guided by this perturbative analysis, we 
propose the following form of the time dependent state in the interaction picture 



|tf (t)> = |* A (t)> \g d ) + \* DS (t)) |0 7 > 



where 



\V A (t)) = C A (t)\A)\0 1 ) + J2 



<?s. + (*) !%,->! + 1 > + ^.-(0lk + )l-i) 



(ni.io) 



(iii.ii) 



and 



\VdsW) = £ [D jt -(t)\ - 1) + D jt+ (t)\ + 1)] \ej) 



(111.12) 



with the initial conditions 



CU(0) = 1; C s± (0) = 0; D ii± (0) = 



(111.13) 



We highlight that I^dsW) describes an entangled state between the spins and the detector. 

The explicit solution for the coefficients with the initial conditions 13|) is provided in 
the appendix. 

The coefficients Dj t ±(t) determine the photodetection probability and display the causal 
nature of the propagation j3l[: the detection time to has to be larger than td = Xd/c, namely 
the time it takes the front of the photon pulse to travel from the NV-center to the position 
of the photodetector. In the experimental setup of ref. [19( the photon travels along a ~ 2 m 
long fiber to the photodetector, therefore td — 7 ns. 

The photodetection probability is obtained as in flIII.8p . and obviously only the state 
\^DA{t)) contributes. The result is a projected reduced density matrix for the spin-qubit 
subpace | ± 1) namely 



p ( i\t) = TTdJ2\^}(4\ xS, Ds(t)}^ D s(t)\ = 



l^- )+ (*)l 2 l + l}(+li + pi,-(*)i 2 i-l}(-ll 



+ (D jt+ (t)D$ t _{t) | + l)(-l|+h.c. 



(111.14) 
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where the coefficients Dj t ±(t) are given in the appendix by (1A.7|) . We now introduce the 
density of states of the photodetector V(u): for any arbitrary function of the detector 
frequencies T{y^) 

/oo 
V{u)F{uj)du ; V(u) = \ K j\K u ~ v 5) • (111.15) 

With the result for Dj t ±(t) given in the appendix (IA.7I) . we introduce 



1 - 


_ gi(w— Q±+i§)(t- 


-t d ) 




tu - £l± + i | 





(111.16) 



in terms of which the projected reduced density matrix at the photodetection time tp in the 
interaction picture becomes 



Pd 



(t D ) = f V(u){\F + (u;t D )\ 2 | + l}(+l| + |.F_l>;f D )| 2 | — X>< — 1 1 

J — OO 

+ 5 P F + (u-t D )r(tu-t D )\ + l)(-l\+h.c.}dtuQ(t D -t d ). (111.17) 

In the narrow width limit T <C fi± the functions J-±{oJ] t) feature sharp peaks at u> = Q± = 
Q ± Aw/2, again we assume that Au <C Q and consequently that T + ~ T_ ~ T/2. We also 
assume a broadband detector whose spectral density is insensitive to the spectral width of the 
emitted photon T and the energy difference between the | ± 1) states Au, namely T>(Q+) ~ 
In particular the correlation function for the broadband photodetector [13, 35] is given 

by 

/oo 
V(u) e 1 ^-^ dcu ~ 2tt V(Q) 5(t - f) . (III. 18) 
-oo 



We can now extract T>(uS) — T>{0) outside the integrals, and using the result ( III. 201) we 
find 



J z a (uj]t D )J z l(uj]t D )dw 



1 _ g— iA ab (tD— to) g— T(t D — tj.) 



1 + i^ 



Going back to the Schroedinger picture at time to and taking r H 
result for the projected reduced density matrix is given by 



A ab = fi a -fi b ;a,b = +, 
(111.19) 

= T_ = T/2 the final 



PD(t 



Dj 



1 - e" 1 * 



0(r) <^| + 1)(+1| + |-1)(-1| 



+ «jr|+l)(-l|e <AwtD 77(r)+h.c.S ; r = t D -t d 



;ill.20) 



where ^(r) is given by ( II 1.27ft with r = tp — t d - 
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Comparing the prefactor of this expression with the total photon number (1II.44|) it is 
clear that the prefactor is just describing the number of photons detected at the retarded 
time to — td and allows the identification of T>(Q) with the detection efficiency. In the 
experimental setup in[l9| this efficiency is <C 1 thus justifying the neglect of the photon 
emission from the decay of the excited states of the detector. The coherence term has a 
simple interpretation: photodetection by filtering the linear polarizations H or V projects 
the spin-qubit-photon entangled state at a time t D into a state similar to that studied 
in section fill A I) effectively disentangling the polarization from the spin degree of freedom 
leaving frequency entanglement only. For Ft ^> 1 the coherence is suppressed by the same 
factor as in the previous case (III.26j) reflecting which path information. 

This result is fully compatible with Glauber's theory of photodetection with an "ideal" 
broadband photodetectorj34-36[, where the detection probability is given by 



P D (t D ) = K / \a ± (t)\E(-\x dl t)E {+ \x d ,t)\a ± {t)) dt = k'^U - e^A 0(r) (111.21) 

1 L J 



here constants [34| and we used eqn. (1A.5|) . Similarly the interference terms are 

given by 



Pi(t D ) = K l t \a + (t)\E^-\x- d ,t)E^\x- d ,t)\a^t))dt= «' 



9(r) 



1 - e- iAuJT e~ Tr 



(111.22) 



These are precisely the terms in the reduced density matrix (1IH.20|) . 

After projection of the photon state into H/V polarization, spin-qubit-photon entangle- 
ment is displayed by projecting on any state of the form 



\M) = -j=[\ + 1) + e^l 



(111.23) 



This is implemented with the reduced density matrix (IIII.20I) by obtaining the conditional 
probability 

PM\ H y(t D ) = Trp D (t D )\M)(M\ . (111.24) 

The non- vanishing coherence in flIII.201) in the basis | ± 1) leads to oscillatory behavior of 
Pm\h,v{^d) a s a function of to- For the state ( 1III.23I) with = and an H projection we 
find for r = to — t d > 



Pm\h{t) 



1 + Re (e iAutD r](r) 



tD — t 



(111.25) 



Fig. 03]) displays the probability ( 1III.25P as a function of r = t D — t d for the experimental 
values reported in ref.0: Aw = 2vr x 122 MHz ; l/Y = 12 ns ; t d = 7 ns. 

This figure reveals the effect of which-path suppression of the coherence: the asymptotic 
behavior of the probability is 



p M \h{t > i/r) 



1 r 



r 



(111.26) 
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x (ns) 



FIG. 3: The probability flIII.25j> for t d = 7 ns , Au = 2vr x 122 MHz ; 1/r = 12 ns. 



Measurement in the //"/V basis results in a post-measurement density matrix that features 
coherence in the qubit basis | ± 1} suppressed by which-path information. This coherence 
was not manifest in the pre-measurement density matrix because of the orthogonality of the 
circularly polarized photon wave packets. 

The reduced density matrix flIII.20p is similar to (III.26p . normalizing so that Pd{ t ) — 
po(r)/Trpo(r) it can be diagonalized in a new basis that differs from (jII47|II.48p by the 
phases multiplying | — 1) and with eigenvalues 



1 ± \ V (r) 



(111.27) 



respectively, leading to the post-photodetection Von-Neumann entropy of entanglement 

S d (t) = -Trp D (r) In p D {r) = -e + (r) lne+(r) - e_(r) lne_(r) . (111.28) 

This post-measurement entanglement entropy is given by Sf (oo) in eqn. (111.5ip asymptot- 
ically for Yt ^> 1 . 



A. Implementing a "Quantum eraser": 

The factor 1 / (l+iAw /V) in the results (1III.20IIIII.22IIII.27P reflects which-path information 
because it suppresses coherence when Au 3> T. It is noteworthy that this suppression 
remains in the final expressions even in an "ideal" broadband photodetector a Id Glauber 
which is insensitive to the photon frequency and with a photodetection correlation function 
oc S(t — t') as discussed above. 

In the experiment in ref-H Au = 2ir x 122 MHz ; F ~ 1/12 ns so that Au/T ~ 
9.2 and there is a strong suppression of coherence because of which-path information 



1/ a/1 + Acu 2 /r 2 ~ 0.11. In this experiment photodetection is carried out with a photode- 
tector with time resolution St ~ 300 ps 1/Aoj to implement a "quantum eraser" 
to "erase" which-path information by introducing an energy uncertainty ~ 1/St ^> Au. 

A simple model for such photodetector can be implemented by modifying the interaction 
Hamiltonian between the detector and the radiation field He>r (lIILlft introducing a " shutter 
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function" S(t) with explicit time dependence, namely 

H DR (t) = Y,[d j -E i P + \x d ;t)\e'l)(g d \e^ t + h.c]s(t) ; P = H;V , (111.29) 

3 

where the only restrictions on the shutter function S(i) are 

m = ln\t D ~ 6t - t - tD (IIL30) 
otherwise 

with the shutter interval 5t such that 

T5t < Au5t < 1 . (111.31) 

This function effectively describes a shutter with a time resolution and amounts to "slic- 
ing" or time-binning the photon wavefunction upon detection. 

A similar procedure of "chopping" the wave function in short time intervals has also been 



advocated as a quantum eraser in ref.|24|. In ref.|30[ a phenomenological damping term 
is added to the right hand side of the equivalent of equations (IA.4|) in this reference, with 
the argument that such damping term describes the coupling of the (single) excited state of 
the detector atom to some reservoir. A "quantum eraser" is implemented in this approach 
by taking the damping constant 7 ^> Aou. While this phenomenological approach seems 
sensible, we consider instead the model of the photodetector with the shutter function S(i) 
introduced above implemented within an ideal broadband photodetector as follows. 
The solution for the coefficients Dj } ±(t£>) are now given by 

D h± (t D ) = -1%. ■ f\%\E { P + \x dl t)\a T (t)) S(t) e^dt, (111.32) 
V 2 Jo 



and the reduced density matrix elements in (1III.14j) become 

Y^D jta {t D )Dl h {t D ) = I f°dt r^t^S^SCO^^I^H^t)^)^^)!^^)) 
• ' z Jo Jo 

/oo 
P(w)e iw( '-'' ) duo ; a, 6 = +,-, (HI.33) 
■00 

where we have used that |cr±(£)) are one-photon wavepackets and only the vacuum con- 
tributes to the intermediate state in the correlation function of the electric field. The last 
term in ( 1III.33I) is the photodetector correlation function (13, 35[ which for a broadband 
photodetector is given by eqn. 18j) . leading to 

Y,D j>a {t D )ir jb {t D ) = 2tt^ f° dt$\t) (omE ( f\x d ,t)E^\x d ,t)\a a (t)) 

~ 27r'^(a b (t D )\E { P \x- d ,t D )E i P + \x-' d ,t D )\a a (t D )) 5t . (111.34) 



where we have used the condition flIII.3ip so that the integrand is constant in the interval 



t D — 5t < t < t D and vanishes outside it. Using the result (IA.5j) we obtain the reduced 
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density matrix in the Schroedinger picture 



Poit 



(r<5t) e- Tr 0(r) 1 1 + 1)(+1| + | - 1)(-1| +5 P (| + 1)(-1| e iAw * d +h.c.) J 



Remarkably, this density matrix describes a pure state, namely 



p D (t D )=M(T) 



+ l) + <Te 



P AQ.-t d \ 



l))(e-' c+td (+l| + ^e 



P -iV.-t dl 



-1 



with the normalization 

N{t) 



(T6tj e- Tr 0(r) ; T = t D -t d . 



(111.35) 
(111.36) 

(111.37) 



It is noteworthy that the quantum eraser has purified the post-measurement reduced 
density matrix. This analysis confirms the experimental results in refs. 19j, |20( and bolsters 
the arguments presented in ref . (20| . 

In the experiment in ref.[l9| after detection the spin-qubit evolves freely in time from tp 
until a time t so that 

p D (t) = M{t) (e in +^\ + 1) + 5 p e in - (t - T) \ - 1)) (e^+ ( '- r) (+l| + 5 P e~ in - {t - T) 

(111.38) 

at which time two microwave pulses resonant with the levels |±) are turned on and transfer 
coherently the state 

' X (111.39) 



\M(t)) = -pie* 1 **! + 1) + e in -*e^| - 1) 
with a fixed phase to the ground state |0), as depicted in fig. (0J. Now we find the total 




FIG. 4: Coherent transfer of the state \M(t)) to the ground state |0) see ref. 19] 
(joint) probability 

A/"(r N 



Pm\h,v{t) = Tr 



p D {t)\M{t))(M(t)\ 



1 ± cos air ) 



a(r) =Awr + f (111.40) 



This result agrees with the joint probability quoted and experimentally confirmed in 
ref. Ji| up to the overall normalization factor and the retardation in the detection time 

T = t D - t d . 
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IV. SUMMARY AND CONCLUSIONS: 

In this article we have studied the dynamics of frequency and polarization entangle- 
ment between photons and a spin-qubit from spontaneous decay in a typical A system with 
non-degenerate lower levels. We addressed in detail how which path information affects co- 
herence, obtained the entanglement entropy for the reduced spin-qubit with frequency and 
polarization entanglement and provided a unified description of the process of spontaneous 
emission and broadband photodetection that is fully causal and allows to include a quantum 
eraser in a consistent manner. 

The main results are the following: beginning with the case in which photon spin-qubit 
entanglement does not involve polarization but only frequency, the reduced qubit density 
matrix obtained from tracing out the radiation bath features oscillatory coherence terms (in 
the qubit basis) that are suppressed by which path information by a factor 1 / a/1 + Aw 2 /r 2 
where Acu is the Zeeman splitting between the lower spin states and T is the linewidth of the 
excited state. In the case in which the spin degree of freedom is entangled with circularly 
polarized photons, the reduced density matrix is a statistical mixture as a consequence of 
the orthogonality of the polarization of the photon states. We obtain the entanglement 
Von-Neumann entropy in both cases and analyze their long time asymptotic behavior. In 
the case in which the spontaneous decay rate is the same to the two lower levels, we find 
that Sf p (t) > Sf (t) where Sf p (t) (Sf (t)) is the entanglement entropy for frequency and 
polarization (frequency only). Focusing on broadband photodetection in the case of fre- 
quency and polarization entanglement, we find that with an ideal photodetector that filters 
photons with horizontal (H) or vertical (V) directions the post-measurement density ma- 
trix describes a mixed state with non-vanishing coherences in the qubit basis. Despite the 
broadband nature of the photodetector described by correlation function oc 5(t — t'), the 
coherences display oscillatory behavior suppressed by which path information just as the 
pre- measurement density matrix in the case of frequency entanglement. 

A "quantum eraser" is implemented within the Glauber model of broadband photodetec- 
tion by including a "shutter function" that effectively time-bins photodetection with a time 
resolution 5t so that T5t <C Au5t <C 1 thereby introducing enough energy uncertainty to 
average out frequency information. We find that photodetection with this "quantum eraser" 
purifies the post-measurement reduced density matrix to a pure state. The resulting joint 
probability for H/V photodetection with projection onto a a superposition of qubit states 



\M(t)) = -^[e in+t \ + l) + e lQ - t e l( t > \ — 1)J is given by ( lIII.40p and agrees with the experimental 
results found in ref . Il9|l . 



Several aspects of the results obtained in this article suggest possible experimental av- 
enues: 1) the dependence on the delay time t d = Xd/c with Xd the position of the photode- 
tector, suggests the possibility of using several photodetectors in coincidence, for example to 
study interference effects, Hanbury- Brown- Twiss correlations or as a complementary vari- 
able to explore coherence as a function of this delay distance, 2) rather than implementing 
a "quantum eraser" with time-binned photodetection, continuous photodetection should 
instead produce a joint probability given by 25|) which displays steps in the coherent 
oscillations (see fig. ©), 3) instead of a "quantum eraser" with time resolution 5t <C 1/Au 
one could consider a "quantum blurrer" with a varying shutter time resolution. This serves 
as a window to let more which path information thereby suppressing the coherence in a 
controlled manner. 

The experimental relevance of the questions studied in this article merit further study 
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perhaps including alternative methods such as those of quantum open systems in terms 



of a master equation [38|, |39J or "quantum jumps" followed by density matrix resetting as 



advocated in ref. 40 



Entanglement and quantum correlations are becoming very important in many ti mely 



aspects of particle physics: in neutrino oscillations [41l. |42j and in CP and T violation |43|, |44 
Recently the entanglement of neutral B-meson pairs produced from the (spontaneous) decay 
of a T(4S) resonance has been exploited experimentally to unambiguously show time- reversal 
violation|45, 461 by tagging individual members of the correlated pairs. Therefore the interest 
on the dynamics of entanglement, the emergence of spontaneous coherence and quantum 
correlations is transcending disciplines and clearly merits deeper understanding. 
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Appendix A: Solutions for the coefficients in eqn. (]III.lipiII.12|) 

The equations of motion for the coefficients in eqn. (1III.llllIII.12j) are obtained from the 
Schroedinger equation in the interaction picture d\^(t))/dt = —iHj(t)\^/(t)) projecting on 
the corresponding states. 

These simplify substantially from the following properties: H A r is the identity in the 
detector space {\e d ), \g d )} and H DR is the identity in the NV-center basis {|^4), | ± 1)}. 

The equations of motion for the coefficients CV ± (t) feature contributions of the form 

(l %± -Tl\g d \H DR \Tl\e d ^) D hT {t) 

arising from the term ^jd* ■ E^(xd,t)\g d )(ej\ in H DR {t). Such term describes the de- 
excitation of the photodetector by spontaneous emission from an excited state \e d ) in which 
the NV-center states | ± 1) are passive, this term is of higher order in dipolar couplings and 
under the assumption of very small detection efficiency as is the case experimentally (see 
below) it will be neglected 3 , leading to the final form of the equations of motion 

iC A {t) = (^O.I^WIl^s-^^J^ + ^O.I^^Il^s+l)^^) (A.l) 
iCfc j+ (t) = {l % ^+l\H AR {t)\A-%)C A {t) (A.2) 
iC M; _(t) = (l£ j+ ;-l|iW*)lA0 7 >CU(t) (A.3) 

Dj,±(t) = • (0,\E P + \x d ,t)\a T (t)} e»*, (A.4) 

3 If necessary, this contribution can be obtained from the unitarity condition (t)\^f (t)) = 1. 
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where the states |<7 T (i)) are given by (III. 411) with (III. 361) evaluated at xq = 0. The solutions 
to eqns. (IA.l|A.2|A.3j) are the same as (III.35IIII.36I) . Upon inserting these solutions in the 
matrix element (IA.4j) . we obtain in the Wigner-Weisskopf approximation 



(1; 



-i^(0 y \E£\x d ,t)\a^{t)) e 



= Kj JE± s? e ^ e *te-n±)(*-*«0 e -§(*-*d) Q(t-t d ) ; t d = — 
V 2ir 



(A.5) C 

where the constants Kj are proportional to dj/xd with proportionality coefficients that result 
from angular and contour integration 4 and 



, P i 1 for P = H 



=F 1 for P = V " 



(A.6) 



From this result we obtain 



Dj,±(t) = iKj 




r± zP Avit. 



S r iv jtd 

2tt t 



1 - 


. gifo— fi±+i|)(i- 











(A.7) 
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